Quantum interference and spin-charge separation 
in a disordered Luttinger liquid 
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We study the influence of spin on the quantum interference of interacting electrons in a single- 
channel disordered quantum wire within the framework of the Luttinger liquid (LL) model. The 
nature of the electron interference in a spinful LL is particularly nontrivial because the elemen- 
tary bosonic excitations that carry charge and spin propagate with different velocities. We extend 
the functional bosonization approach to treat the fermionic and bosonic degrees of freedom in a 
disordered spinful LL on an equal footing. We analyze the effect of spin-charge separation at fi- 
nite temperature both on the spectral properties of single-particle fermionic excitations and on 
the conductivity of a disordered quantum wire. We demonstrate that the notion of weak localiza- 
tion, related to the interference of multiple-scattered electron waves and their decoherence due to 
electron-electron scattering, remains applicable to the spin-charge separated system. The relevant 
dephasing length, governed by the interplay of electron-electron interaction and spin-charge sepa- 
ration, is found to be parametrically shorter than in a spinless LL. We calculate both the quantum 
(weak localization) and classical (memory effect) corrections to the conductivity of a disordered 
spinful LL. The classical correction is shown to dominate in the limit of high temperature. 

PACS numbers: 71.10.Pm, 73.21.-b, 73.63.-b, 73.20.Jc 



I. INTRODUCTION 



Interacting electrons in one dimension (ID) are 
a paradigmatic example of the strongly correlated 
fermionic systems. Electron-electron (e-e) interactions 
drive the ID system into a non-Fermi liquid state known 
as the Luttinger liquid (LL) (for review see, e.g., Refs. [0, 
In recent years, progress in nanofabrication 
technologies has made it possible to manufacture a vari- 
ety of single- and few-channel quantum wires connected 
to the electric leads and to perform systematic trans- 
port measurements on the very narrow wires. The lat- 
ter include single-wall carbon nanotubes^ ' 10 ! 11 ' 12 ! 13 ' 14 , 
semiconductor-base d 15 ! 16 ' 17 : 18 and metallic^ quantum 
wires, polymer nanofibers 2 ^, as well as quantum Hall 
edge states! 21 ! 22 The LL nature of these strongly corre- 
lated quantum wires has been supported by a wealth of 
experimental findings^a. n i 12 ! 13 ! 14 . 15 i 16 ! 19 i 20 ' 21 i 22 An- 
other class of strongly-correlated quantum wires that has 
recently attracted a lot of interest is ultracold atomic 
gases confined to ID geometry, for review see Ref. |23|. 

Mesoscopic physics of strongly correlated electrons is 
one of the most important and promising directions of 
current research on ID electron systems. Recent trans- 
port measurements on carbon nanotubes 2 ^ reported both 
sample-dependent conductance fluctuations and strong 
magnetoconductivity, in qualitative similarity to the 
mesoscopic phenomena in higher-dimensional disordered 
electron systems. In Ref. [2J, the sample size ~ 1 /im 
was of the order of the mean free path limited by im- 
purity scattering. Electron transport through the nan- 
otubes displayed, therefore, features characteristic of 



the crossover from ballistic conduction to a disorder- 
dominated regime. On the other hand, in the past few 
years techniques to grow nanotubes of size up to the mil- 
limeter scal o 25 ! 26 — much larger than the typical value of 
the disorder-induced mean free path — have been devel- 
oped. First transport measurement a 25 ' 26 on the ultra- 
long nanotubes provided evidence for disorder-induced 
diffusive motion of electrons in a wide range of tempera- 
ture. Altogether, these advances have paved the way for 
systematic experimental study of interference-induced lo- 
calization phenomena in ID. 

The theory of weak localization (WL) in a disordered 
LL of spinless electrons was developed in Refs. and 
l28l (for recent advances in the ballistic cr-model frame- 
work, see also Ref. l29t ). In the limit of weak interac- 
tion between electrons, the phase breaking length 1$ that 
sets up the infrared cutoff of WL was shown to obey 
1$ ~ a -1 ^/) 1 / 2 , where a <C 1 is the dimensionless in- 
teraction constant, It ~ v/T the thermal length, v the 
Fermi velocity, T the temperature, I the elastic mean free 
path. At sufficiently high temperatures, when 1$ <C /, 
the system is in the WL regime.— The WL correction 
Actwl to the Drude conductivity ctd behaves then as 



-(V0 2 ln(//W 



\l T /l)Ha 2 l/h 



Note that the WL dephasing length of spinless electrons 
originates from the interplay of interaction and disorder: 
both ingredients are necessary to establish a nonzero de- 
phasing rate. As a result, 1$ in the WL regime is much 
longer than the total decay length of fermionic excitations 



a It- The latter is the dephasing length relevant 



to the damping of Aharonov-Bohm oscillation o 27 ! 28 ' 31 ! 32 
and to the smearing of a zero-bias anomaly in the tun- 
nelling density of states^ 3 - 



2 



In a spinful LL, the most prominent spin-related man- 
ifestation of non- Fermi liquid physics is a phenomenon of 
spin-charge separation (SCS) (see, e.g., Refs. [3(301300, 
7). The essence of the SCS in the LL model is that the 
spin and charge sectors of the theory in a bosonic repre- 
sentation are completely decoupled from each other and 
characterized by different interaction coupling constants. 
Correspondingly, the elementary bosonic excitations car- 
rying spin and charge propagate with different veloci- 
ties, independently of each other. Experimentally, the 
effect of the SCS on the spectral properties of a LL (as 
measured in electron tunnelling experiments) has been 
studied in Refs. H5IIT7I . In the last few years, much at- 
tention has been given to the so-called "spin-incoherent 
regime" 34 in ID systems with strongly different spin and 
charge velocities i 35 i 36 ' 37 For recent transport measure- 
ments in ballistic quantum wires which show signatures 
of the spin-incoherent behavior, see Ref. [IK and refer- 
ences therein. 

While the difference in the spectral properties of spin 
and charge collective modes is not at all specific to ID 
and is also characteristic of higher-dimensional Fermi 
liquids f 3 ^ the peculiarity of ID is that the "factoriza- 
tion" of the bosonic modes modifies the single-particle 
fermionic properties in an essential way— Our purpose 
is to investigate how the SCS affects the quantum inter- 
ference phenomena in a disordered LL. In particular, we 
analyze the dynamical properties of fermionic excitations 
at finite temperature and employ the results of this anal- 
ysis to calculate the WL correction to the conductivity 
of a spinful LL. The problem is rather nontrivial con- 
ceptually since the spin and charge degrees of freedom 
that constitute the electron in a LL acquire different ve- 
locities and the very notion of a specific quasiclassical 
electron trajectory characterized by a certain velocity — 
conventionally invoked in a description of WL — becomes 
ambiguous. 

We show below that, despite the intricate nature of 
a single-electron motion in a spinful LL, the basic no- 
tions of WL remain applicable even when the SCS is 
incorporated in the calculation. However, the spin de- 
gree of freedom has dramatic consequences for the effects 
of e-e scattering. Most importantly, the decay rate of 
fermionic excitations is strongly enhanced in the presence 
of spin: the single-particle length l ee ~ a -1 It becomes 
parametrically shorter than for spinless (spin-polarized) 
electrons. Moreover, in contrast to the spinless case, the 
WL dephasing length 1$ ~ a' 1 It becomes of the or- 
der of lee- As a result, also the WL dephasing turns 
out to be much stronger than without the SCS being 
included. The temperature at which spinful electrons 
get strongly localized (It /I ~ ct) is therefore much lower 
(for small a) than for spinless electrons. Furthermore, 
we find that the WL correction to the conductivity is 
given by Ad WL /(T fl ~ -(/^/Z) 2 ~ -~a~ 2 (l T /l) 2 , showing 
a much faster temperature dependence of Acjwl than in 
the spinless case. 

We also demonstrate that a classical "memory effect" 



(ME) in the electron scattering off disorder contributes to 
the T dependence of the conductivity. Moreover, it gives 
the leading (larger than Acjwl) correction to the Drude 
conductivity in the limit of high T. The obtained ME 
contribution Actme/c_d ~ —It /I is essentially not related 
to e-e interactions and exceeds Actwl only when the lat- 
ter is sufficiently suppressed by the interaction-induced 
dephasing. Specifically, |Acjme| 3> |Actwl| for It /I <C a 2 
(i.e., for l^/l < a). Otherwise, |Acr WL | > |Actme|- 
Technically, the ME manifests itself in the same set of 
Feynman diagrams for the conductivity as the WL, so 
that we actually treat the two effects — the essentially 
classical ME and the essentially quantum WL — on an 
equal footing. What distinguishes them from each other 
is that the main contributions to Actwl and Acme come 
from scattering on different impurity configurations. The 
WL correction stems from scattering on rare, compact 
three-impurity complexes in which the characteristic dis- 
tance between all three impurities is of the order of the 
single-particle length l ee . The ME correction is associ- 
ated with impurity configurations in which two of impu- 
rities are located very close to each other, with a charac- 
teristic distance between them of the order of the thermal 
length It- 

The remainder of the paper is organized as follows. In 
Sec.|TT]we formulate the model of a disordered LL, discuss 
the basic physics of the SCS relevant to the e-e scattering, 
and describe the method of functional bosonization used 
in our calculation. Scction HTTl is devoted to an analysis of 
the spectral properties of fermionic excitations in various 
representations; in particular, in a "space-energy repre- 
sentation" employed for the calculation of the conductiv- 
ity. In Sees. IIVI and fVl we evaluate the WL and ME cor- 
rections to the conductivity, respectively, by means of the 
functional bosonization method. In Sec. IVII we present a 
complementary analysis based on the more conventional 
path-integral approach. Our results are summarized in 
Sec. EHl 



II. MODEL AND METHOD 

We begin by formulating the model of a disordered 
spinful LL in Sec. Ill Al In Sec. Ill B I we present a simple 
argument which demonstrates the peculiarity of ID ge- 
ometry in that the spin degree of freedom affects the rate 
of e-e scattering in the LL in a crucial way. Section III CI 
is devoted to an overview of the functional bosonization 
method. 



A. Disordered Luttinger liquid 

Throughout the paper we consider a single-channel in- 
finite quantum wire. Linearizing the dispersion relation 
of electrons about two Fermi points at the wavevectors 
k = ±kp with the velocity v, the Hamiltonian of a clean 
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LL is written as (h = 1) 

kfJ,(T 

/Acer' 

(2.1) 

Here ^^(fc) are the electron operators at the wavevector 
k, the index fi = ± denotes two branches of chiral exci- 
tations (right and left movers), and a =|, J, stands for 
two spin projections. The e-e interaction enters Eq. (|2.1[) 
through the coupling constants g± and gi- These describe 
forward e-e scattering with small momentum transfer 
(much smaller than hp) between electrons from the same 
(34) or different (32) chiral branches. We assume that the 
e-e interaction is short-ranged and represent the interac- 
tion part of the Hamiltonian in terms of the local in space 
electron density operators n fia (x) = ip^ix)^^^). The 
parameters of the Hamiltonian with the linearized disper- 
sion should be understood as effective (phenomenologi- 
cal) couplings of the low-energy theory, which include 
possible high-energy renormalization effects, similar to 
Fermi-liquid theory. 

The LL Hamiltonian (|2.ip does not contain the term 

Hhs = \ J dx ^\o^-w » ( 2 - 2 ) 

which describes backward e-e scattering with large mo- 
mentum transfer resulting in a change of chirality fi. If 
one begins with a microscopic model of electrons inter- 
acting via a finite-range externally screened Coulomb po- 
tential, the constants (72.4 in Eq. (|2.ip and g\ in Eq. (|2.2j) 
are related to the Fourier transforms of this potential 
at zero and 2kp momenta, respectively. The forward 
scattering dominates provided that the radius of external 
screening d (e.g., the distance to a metallic gate) is much 
larger than kp . We assume that this is the case and 
neglect f/bs throughout the paper below 4^ Treating the 
Coulomb potential in Eq. (|2.1[) as short-ranged is legit- 
imate for scattering processes with momentum transfer 
much smaller than d . This same scale <i _1 fixes the 
ultraviolet momentum cutoff in our formulation of the 
low-energy theory. 

The only source of electron backscattering in our model 
is thus a static random potential U (x) due to the presence 
of impurities. We assume that fluctuations of U (x) are 
Gaussian and characterized by the correlation function 
(U(x)U{x')) = 5{x - x')v 2 /2l ("white noise"). Here l 
is the transport elastic mean free path in the absence 
of interaction. The disorder is considered to be weak, 
kplo 3> 1. The disorder-induced backscattering term in 
the Hamiltonian is given by 

H imp = ^2 [dx (U; V^>- + U b ^I a lM , (2.3) 



where the backscattering amplitudes U b (x) are correlated 
as (U b (x)U£(x')) = (U(x)U(x')) and (U b (x)U b (x')) = 0. 
Forward scattering off impurities can be gauged out in 
the calculation of the conductivity^^ and will therefore 
be neglected from the very beginning. The total Hamil- 
tonian H that defines our model of a disordered spinful 
LL is thus 

H = H Lh + H imp . (2.4) 

Throughout the paper we consider a quantum wire 
with spin and chiral channels not separated spatially in 
the transverse direction, so that the constants 92 = .94 = 
g are spin-independent and equal to each other. The 
plasmon velocity u for the spinful case then reads 

u = v/K p = v{l + 2g/nv) 1/2 , (2.5) 

with K p being the Luttinger constant in the charge sec- 
tor, whereas the velocity of elementary spin excitations is 
equal to v. It is convenient to characterize the strength of 
e-e interaction by the dimensionless coupling constant^ 
a = (1 - K 2 p )/(l + 3K 2 ), which in the limit of weak 
interaction a <^ 1 is written as 

a~ (l-K p )/2 ~g/2-KV . (2.6) 

B. Why spin matters 

To qualitatively understand the nature of dephasing of 
fermionic excitations in a spinful LL, it is instructive to 
recall the perturbative expansion&22 of the self-energy of 
the single-particle Green's function in the limit of weak 
interaction a -C 1 and discuss the e-e scattering rate at 
the Goldcn-rulc level, first in the absence of disorder. For 
the spinless case, such an analysis has been performed 
in Refs. (23111, and H— see also Refs. H^lllEillll and 
Refs. l27l l 32l for a closely related calculation of the tempo- 
ral decay of the single-particle Green's function for spin- 
less and spinful electrons, respectively. Since the physics 
of dephasing is governed by inelastic e-e scattering, a nat- 
ural first step is to calculate at lowest (second) order in 
a the e-e scattering rate t^} given by the imaginary part 
of the self-energy. 

The Golden-rule expression for the e-e collision rate 
reads 

— = f dude' K(lo) (ftJe>ft +a + h-uftU+J) , 

^ee vW J 

(2.7) 

where 

K(u)=r) s [K% + [u)+K*_{u)} + K F {uj) (2.8) 

is the kernel of the e-e collision integral and f e is the 
Fermi distribution function, f^ 1 = 1 — f e . In Eq. (|2.8[) . 
the Hartree terms K+ + cx g\ and K?_ cx g\ are related 
to scattering of two electrons from the same (++) or dif- 
ferent (-1 — ) chiral branches, respectively, K F = —K¥, 
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is the exchange counterpart of K¥,, r/ s is the spin de- 
generacy, rj s — 1 for the spinless case and rj s = 2 for the 
spinful case. 

To order 0(a 2 ), the Golden-rule scattering rate and 
the self-energy on the mass shell coincide with each other. 
At the Fermi level (e = 0) we have 

= -2 [ V s(lm^ + + Im£«_) + ImS F ] , (2.9) 

where the Hartree terms are given by 



ImE 



H 



-± 



= a v diouj ( coth • 



2 J V 2T 

x / dq5(u) — vq)5(u> =F vq) 



tanh • 



2T") 
(2.10) 



and the exchange term T, F = — Peculiar to ID are 
highly singular contributions to K(u) related to scat- 
tering of electrons moving in the same direction. One 
sees that the contribution of £+ + contains a 5-function 
squared and thus divergesi 6 ' 27 i 28 The divergency of the 
perturbative expression for the probability of scattering 
of two electrons of the same chirality simply means that 
the energy and momentum conservation laws for this kind 
of scattering give a single equation lu — vq = 0. 

For spinless (spin-polarized) electrons, the divergency 
in ImE^ , is canceled by the same divergency in the ex- 
change term. The remaining term ImE^_ yields 27 - ' 28 i 32 



1 



r GR 



(0) 



a 2 nT 



(2.11) 



Note that for scattering of electrons from different chi- 
ral branches on each other, the energy and momentum 
conservation laws lead to two equalities: ui ~ vq = and 
uj+vq — 0, which combine to give to, q = for allowed en- 
ergy and momentum transfers. This "quasi-elastic" 27 i 28 
character of e-e scattering is a peculiarity of ID: in higher 
dimensionalities, the characteristic energy transfer that 
determines t~ in a clean system is of order T. 

For spinful electrons, the Fock contribution cancels 
only the part of the Hartree term ImE^ + that comes 
from interaction between electrons with the same spin. 
The divergent second-order Hartree term that arises from 
interaction between electrons with opposite spins remains 
uncompensated. This indicates that the main contribu- 
tion to r~ is now related to scattering of electrons from 
the same chiral branch. Thus, already the perturbative 
expansion demonstrate a 6 ' 27 > 28 a qualitative difference be- 
tween the cases of spinless and spinful electrons. 

In fact, for spinful electrons, the perturbative expan- 
sion of T^ 1 in powers of a is diverging in the clean limit 
at each order. We will analyze the finite-T damping of 
the single-particle Green's function for a <C 1 in Sec. IIIII 
Here, we stick to the calculation of within a "gen- 

eralized Golden-rule" scheme. The term "generalized" 
means that we go beyond second order in a by introduc- 
ing the dynamically screened e-e interaction V{uj, q) — 
which is exactly^ given by the random phase approxima- 
tion (RPA) [see Eq. (|2.25p below] . The second (5-function 



in the integrand of Eq. (|2 . 10[) comes precisely from the 
imaginary part of the retarded propagator ImV(q,uj) if 
one takes the propagator at second order in a, 



lmV(q 7 uj) ~ — (2na) 2 UJv[ 5(u> — vq) + S(u + vq) 



(2.12) 



Using the full RPA propagator ImV (q, ui) cx [ S(uj — uq) + 
5(lu + uq) ], the leading at a < 1 Golden-rule expression 
for the e-e scattering rate of spinful electrons is written 
as 



1 



r GR 



(0) 



2na 2 vT duj dq5(ui — vq)S(u>~uq) , (2.13) 



which only differs from -2Im£^ + in Eq. (|2~TU|) in that 
one of the (^-functions has a shifted velocity v — * u. Note 
that, similarly to the spinless case, l/r^ R (0) in Eq. (|2.13p 
is determined by u, q = 0. The relative shift between the 
arguments of the <5-functions makes the expression for 
l/r£ R (0) finite: 



1 



r GR 



(0) 



2na 2 T ■ 



\a\irT , 



(2.14) 



where we used Eq. (|2.5|) for u — v ~ 2av for small a. 
Remarkably, the e-e scattering rate for spinful electrons 
turns out to be of first order in a, in contrast to the 
spinless case (|2.11l) . where it is of order a 2 . 

As we will see below in a more consistent treatment 
which does not rely on the generalized Golden-rule ap- 
proach, the scattering rate l/r ( £ R (0) gives a characteris- 
tic decay rate for single-particle excitations and also the 
characteristic dephasing rate for WL. What the above 
consideration teaches us is that in ID the spin degree 
of freedom strongly enhances the e-e scattering rate for 
weakly interacting electrons. The parametric difference 
between the spinless and spinful cases is in stark con- 
trast to higher dimensionalities, where taking spin into 
account typically yields for relaxation rates in a clean 
system only numerical factors of order unity. 



C. Functional bosonization 

The method we use here to study the quantum interfer- 
ence in a disordered spinful LL is functional bosonization. 
It was introduced for the clean LL model in Refs. I45H461 
and further developed in Refs. |47||48||49||50| . In the ear- 
lier wor k 27 i 28 by three of us, the functional bosonization 
framework was extended to deal with disordered prob- 
lems and applied to study the transport properties of 
a disordered spinless LL. In this subsection we present 
a brief outline of the formalism (for more details see 
Sec. VII in Ref . 28) , highlighting the differences between 
the spinless and spinful cases. 

In contrast to "full bosonization" , conventionally used 
for a theoretical description of the LL, the functional 
bosonization technique preserves both fermionic (elec- 
trons) and bosonic (collective excitations - plasmons, 
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spinons) degrees of freedom. This feature of the 
method is of great advantage when one has to deal 
with interacting problems which are most naturally de- 
scribed in terms of fermionic excitations, e.g., quan- 
tum interference (Refs. [27ll28l and the present work) or 
nonequlibriu m 33 ' 5 1 > 52 phenomena in a LL. In particu- 
lar, the functional bosonization allows for a straightfor- 
ward treatment of e-e interaction while residing in the 
fermionic basis, which is especially cost-efficient in the 
disordered case. 

The key steps in setting up the formalism for a LL at 
thermal equilibrium are as follows: 

• a conventional Hubbard-Stratonovich decoupling of 
the four-fermion interaction term in the Matsub- 
ara action is performed by means of introducing a 
bosonic field tp(x, r); 

• interaction of fermions with the field tp is gauged 
out by means of a local transformation (/i = ±) 

4>p,o{x,t) -> tp fl(7 (x,T)exp[i6 IJ .(x,T)], (2.15) 

where the phase 9^{x,t) obeys 

(d T -i(ivd x ) 6^{x,t) = ip(x, t). (2.16) 

This transformation completely eliminates the cou- 
pling between the fermionic and bosonic fields from 
the action (this property is peculiar to ID); 

• upon this transformation, the bosonic part of the 
action remains Gaussian. It is this point at which 
the peculiarity of the LL model — the exactness of 
the RPA — comes into play. The correlation func- 
tion of the field tp is given by 



(<p(x,t)<p(0,0))=V(x,t), 



(2.17) 



where V(x, r) is the dynamically screened interac- 
tion [see Eq. ([2~25]) below]; 

• an arbitrary time-ordered fermionic average is ex- 
pressed through a product of free electron Green's 
functions and Gaussian averages of the phase fac- 
tors cxp[i6^(x, t)] (taken at different space-time 
points). The bosonic averages are represented in 
terms of the correlation functions 

b^{x, t) = ( [0„(o, o) - e^x, t)} e u (o, o) ) 

(2.18) 

which are related to the Fourier component 
V(q,i£l n ) of the interaction propagator (|2.17[) as 

B +± (x.r) =T^f^ 

n 

(vq-iQ n )(±vq-ift n ) 
B — (x,t) = B ++ (-x,t) , 

B_+(x,r) = B+_(x,t) . (2.19) 



Here Q r , 
quency; 



2imT is the bosonic Matsubara fre- 



• while calculating observables (closed fermionic 
loops), e-e interaction is completely accounted 
for by attaching the fluctuating gauge factors to 
backscattering vertices. If the number of fermionic 
loops is larger than one, each of them has to contain 
at least one pair of backscattering vertices in or- 
der not to be disconnected (Wick's theorem for the 
functional bosonization diagrammatic technique). 



-i6 + (l) 



+i6 + (2) 



.-B+H-fe-Xi-Vti) 



+ 




FIG. 1: The Green's function of a right mover propagating 
between space-time points 1 = (xi,n) and 2 = (x2,T2) be- 
fore (a) and after (b) averaging over fluctuations of the gauge 
factors. Solid line: the bare Green's function. The wavy 
lines at the end points represent the factors exp[— i8+(l)] and 
exp[i#+(2)]. The wavy line connecting points 1 and 2 denotes 
averaging over fluctuations of #+(1) and #+(2). 

As a simple example, consider the single-particle 
Green's function for, say, a right mover G f +(a;,r). Upon 
gauge transformation (|2.15[) . one gets the free Green's 
function g+(x, t) dressed by two phase factors as shown 
in Fig. QJa). Pairing of the two bosonic fields yields 
[Fig. Hp))] 



G+(x,t) =g + (x,r) exp [-B ++ (x, t)] 



(2.20) 



More complex quantities are calculated in a similar 
way. Each impurity backscattering vertex at space- 
time point N generates a phase factor of the type 
exp{±i[8 + (N) - 0-(N)]}, as illustrated in Fig. H Upon 
averaging, the phase factors are paired in all possible 
ways (Fig. [3]). In closed fermionic loops, the correlators 
B^ v (x,t) [Eq. (|2.18p ] only appear in the combination 

M(x, t) = B ++ {x,t)+ B—(x, t) - 2B+- (x, r) . 

(2.21) 

As a result, each pair of backscattering vertices at points 
(xn,tn) and {xn',tn') contributes either the factor 



Q{x 7 t) = exp[M(x,T)] 



(2.22) 



where x = xn — xn>, t — tn — tjv', or Q^ 1 (x,t), de- 
pending on whether chirality of incident electrons at the 
vertices is the same (Q) or different (Q _1 ). 

The RPA dynamically screened interaction V(q, i£l n ) 
obeys 



V-ifaiSln) = g- 1 +U(q,in n ) 



(2.23) 
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-ie + (i) 



i[9 + (2)-9.(2)] 
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FIG. 2: Backscattering of a right mover off an impurity (de- 
noted by a cross) at point 2. The impurity vertex is dressed 
by a local gauge factor exp{i[#+(2) — #_(2)]} which contains 
two fluctuating fields of different chirality. 



B + _(3,l) 




B + _(2,l) 



B + _(3,2) 



B + .(2,2) 



FIG. 3: Backscattering off an impurity as shown in Fig. [2] 
after averaging over fluctuating bosonic fields. Each of the 
wavy lines represents a factor of the type exp(±.B M „). 



where iVL n ) is the polarization operator. In a spinful 
LL, the latter is written as 

IL(q,iQ n ) = A . ( 2 - 24 ) 



which gives 



ttv v 2 q 2 + fl^ 



v 2 a 2 + n 2 



(2.25) 



with u from Eq. (|2.5p . 

As will be seen below, it suffices, when calculating the 
dephasing rate for weak localization, to deal with the 
ballistic interaction propagator (|2.25p which does not 
include backscattering of electrons off disorder. This 
should be contrasted with the spinless case, where the 
dephasing of localization effects is absent altogether un- 
less the disorder-induced damping of V(q,ifl n ) is taken 
into account^^ ("dirty RPA"). 

Substituting Eq. ([2~25]) in Eqs. (|2~T9|) yields 



B++(x,t) = -^ln7?(x,r) - ^-lns(x,T) 
B+-(x,t) = — - lns(x,T) , 



(2.26) 



where 
T)(x,t) 



(ttT/A) 2 



sinh[7rT(x/u + it) ] sinh[7rT(a;/u — it) 
v sinh[7rT(x/u + it)] 



u sinh[7rT(x/?i + it)] 



and the constants a b and a r are given by 



a b 



(u — v) 
2uv 



2uv 



(2.27) 



(2.28) 



Inspecting Eqs. (|2.26p and (|2.27p . we see that there is 

an extra factor of 1/2 in front of both B ++ and as 

compared to the spinless^ case. It is this factor that is 
responsible for the SCS. 

The exponents a b and a e — a r + at determine a power- 
law suppression of the tunneling density of states (zero- 
bias anomaly) for tunneling in the bulk and in the end of 
a LL, respectively (see, e.g., Refs. [1.2.3.4.5.6,3). In the 
rest of the paper, we will treat the interaction strength 
a as a small parameter. The hierarchy of the constants 
(|2.28p is then as follows 



^ < a r « 1 



(2.29) 



since a b is quadratic in a, whereas a r is linear. Using 
this hierarchy will greatly simplify the calculation below. 



III. SINGLE-PARTICLE SPECTRAL 
PROPERTIES AND SPIN-CHARGE 
SEPARATION 

In this section, we use the functional bosonization for- 
malism to study the single-particle spectral properties 
of electrons at finite temperature. We begin with the 
space-time representation in Sec. IIII Al Then we trans- 
form to a "mixed" space-energy (Sec. IIII Bp and the 
momentum- energy (Sec. lIIl"Cl) representations, using ap- 
proximations appropriate in the weak-interaction limit 
The analysis of the various representations of the 
single-particle Green's function will serve as a starting 
point for the calculation of the WL and ME terms in the 
conductivity in Sees. IIVI andfVl 



A. Green's function in the (a?, r) representation: 
Weak-interaction approximation 

In the absence of interaction, the single-particle 
Green's function of right (+) and left (— ) movers g±(x, r) 
is given by 



iT 

9±(x,t) = Ttt: 3 



1 



2v sinh[7rT(x/i> ± it)] 



(3.1) 



Plugging Eqs. (j2~^6l) . (|2~2Tp . and (EU into Eq. (|2T2T?j) . we 

have for the Green's function of right movers in a spinful 
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LL: 



G+(x,r) 



2tt^/uv 



1/2 



sinh[7rT(x/u + ir)] sinh[7rT(x/u + ir)] 

ttT/A ttT/A 1 ° b/4 

sinh[7rT(:r/u + it)] sinh[7rT(a;/u — ir)] J 

(3.2) 

in agreement with the result obtained by conventional 
bosonization (see, e.g., Refs. dlMlCIIIG) an d by purely 
fermionic methods (see, e.g., Ref. |l|). The Green's func- 
tion of left movers G-(x,t) = G + (—x,t). The ana- 
lytical structure of G+(x,t) in the complex plane of r 
(0 < Rer < l/T) is shown in Fig. |4] for positive x. There 
are three branch points at r = ix/u, ix/v, and —ix/u. 
One way to choose branch cuts is shown in the top left 
panel of Fig. [4] those starting at r = ix/u and ix/v are 
sent upwards, whereas that starting at r = —ix/u is sent 
downwards. If a <C 1, the first two cuts are much differ- 
ent from the third one. In the limit of small a, the cut 
that connects the points t = ix/u and ix/v corresponds 
to an almost square-root singularity, so that the main 
change the Green's function experiences when crossing 
this cut is a change of sign ( "strong cut" ) . On the other 
hand, the cut that goes from —ix/u to —zoo is "weak" in 
the sense that the discontinuity of G+{x,t) across this 
cut is proportional to aj, ~ a 2 <C 1. Similarly, cross- 
ing the axis of imaginary r between r = ix/v and ioo is 
associated with a weak discontinuity. 

The main approximation we make in this paper con- 
sists in sending ah to zero everywhere in the calculation 
while keeping the effects of leading (linear) order in the 
interaction strength. That is, below we retain the differ- 
ence between u and v, Eq. (|2.5[) . 



(3.3) 



v(l + 2a) , a < 1 



as the only effect of e-e interaction^ The Green's func- 
tion G+(x,t) then reads 



G+(x,t) 



2n^/uv 
ttT 



nT 



smh[?rT(x/v + ir)] sinh[7rT(a;/it + ir)] 



1/2 



(3.4) 

The velocities u and v in Eq. (|3.4p coincide with the ve- 
locities of the elementary collective excitations (plasmons 
and spinons). The appearance of the two velocities in the 
single-particle correlator signifies SCS. Within the ap- 
proximation (|3.4p , the two velocities enter the fermionic 
Green's function in a symmetric way. The analytical 
structure of G+(x, t) in Eq. (|3.4[) is simplified to a single 
square-root cut between the points r = ix/u and ix/v, 
as illustrated in the center top panel of Fig. 0] We will 



i xl v 
i xl u 



-ix/u 



i xl v 
ix/u 



i xl v 



i xl u 







P 



i x/v 



i x/u 







P 



FIG. 4: Left top panel: the analytical structure of the right- 
mover Green's function G+(x,t) in the complex plane of the 
Matsubara time r. The bold solid and the dashed lines rep- 
resent "strong" and "weak" branch cuts, respectively [see the 
text below Eq. (|3.2|l ]. Central top panel: within the approx- 
imation (|3.4p . only the square-root branch cut between the 
points r = ix/u and ix/v survives. Last three panels: the 
contour transformation used to calculate the Green's func- 
tion in the space-energy representation in Sec. IIIIBI 



use the approximation (|3.4p , which captures the essential 
physics of SCS, throughout the paper below. 

It is instructive to compare Eq. (|3.4[) with the Green's 

(x,t) (see, e.g., Refs. [l], 



function of spinless electrons 



tunction ol spir 

2.3.4.5.6.7.28). Since the correlator B ++ ( x,t) in the 
spinless case is twice as large, G+(x,t) is given by 



ttT 



2ttu sinh[7rT(x/u + ir)] 
ttT/A ttT/A 1 ab/2 



sinh[7rT'(x/u + ir)] sinh[7rT(a;/u — ir)] 



(3.5) 



One sees that in the absence of spin, interaction leads to 
a replacement of the velocity v — > u in the bare Green's 
function and generates two brunch cuts: one with an ex- 
ponent close to 1 and the other with the small exponent 
a.h/2. It follows that for spinless electrons the approxi- 
mation at, — * 0, analogous to Eq. (|3.4p , would eliminate 
all dephasing effects — since the latter only originate from 
the factors in the second line of Eq. (|3.5p . By contrast, 
dephasing in the spinful case arises already at order 0(a), 
as will be shown below. Consequently, the approxima- 
tion (|3.4p allows us to obtain — in a controllable way — 
analytical results valid in the limit a <C 1. 

Let us now identify two important spatial scales. For 
this purpose, we perform the Wick rotation t — > i(t + i0) 
inEq. For large \x/u-t\, \x/v-t\ > l/T, Eq. I^g]) 

yields 

G+(x,it) oc exp[-TrT(\x/u-t\ + \x/v -t|)/2] . (3.6) 
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Within the interval x/u < t < x/v, the Green's function 
given by Eq. (|3.6|) decays as 



G+(x, it) oc exp (— x/2l ee ) , 



(3.7) 



independently of t, whereas outside this interval the 
Green's function is suppressed much more strongly; in 
particular, at t = 0: 



G + (x, 0) oc exp (— x/21t) , 



(3.8) 



merely due to the thermal smearing. Here we have intro- 
duced 



and 



2ttT 



2ttT 



(3.9) 



(3.10) 



which are the length scale of spatial decay of fermionic ex- 
citations due to e-e interaction and the "thermal smear- 
ing length" , respectively. The velocities u± are given by 



1 _ 1 / 1 1 

u± 2 \v u 



(3.11) 



The length l ee has also been termed the Aharonov-Bohm 
dephasing length. 27 , i 28 i 31 i 32 Note that this length agrees 
with the Golden-rule estimate ()2.14|) . up to a numerical 
factor. 

For a <C 1 we have 



l T ^ v/2nT, l ee ~ It /a, 



(3.12) 



i.e., l ee for weak interaction is much longer than It- In 
Sees. II Vi and V, when considering the system in the pres- 
ence of disorder, there will appear one more character- 
istic length scale: the electron mean free path due to 
backscattering off impurities I . We will assume that T is 
sufficiently large, so that l ee <C I. As will be seen below, 
this condition means that the disordered system is in the 
WL regime. For lower temperatures, strong localization 
sets in.— Altogether, the hierarchy of length scales in our 
problem is 



It <l ee <l 



(3.13) 



B. Green's function in the (x,e) representation 

We now turn to the single-particle Green's function 
in the space-energy representation, which is obtained by 
Fourier-transforming G+(x, r) with respect to r. Within 
the small-a approximation (|3.4[) . the only singularity of 
G+(x,t) is a branch cut between r = \x\/u and \x\/v in 
the upper or lower half-plane of r depending on the sign 
of x. Since G + (x, r) in this approximation is analytical in 
one of the half-planes of r, its Fourier transform vanishes 
for e„ < if x > or for e n > if x < 0. For both e n 



and x positive, we transform the contour of integration 
as shown in Fig. 2J Closing the contour upwards, the 
integral along the real axis of r is represented as a sum of 
two integrals along the imaginary axis at r — +0 and r — 
1/T — 0. In view of the periodicity of G+(x, r) in r, the 
sum gives the integral along the contour around the cut. 
Closing similarly the contour of integration downward if 
both e n and x are negative, we get 

l/T 

dr exp(ie„T) G + (x, t) = G r + (x,ie n ) - G\{x,i£ n ) . 

(3.14) 

Here e n = 2-k(ji + h)T is the fermionic Matsubara fre- 
quency, 

G r + (x,ie n ) = 6(s n )9(x)<g(x,e n ) , 

G%( Xl ie n ) = 9(-e n )9(-x)&(x,£ n ) , (3.15) 

and the function Sf(x, e n ) depends on the absolute values 
of the coordinate and energy: 

&(x,e n ) = —== exp(-\e n x\/u) / dt 
W uv Jo 

exp(-|e„|t) 

[sinh(TrTi) sinh(27rT|x|/u_ - nTt)} 1 / 2 ' 

(3.16) 

For left movers, G r _! a (x 1 ie n ) = G+ a {— x, ie n )- Integration 
in Eq. (|3.16[) yields (in the rest of the subsection, let both 
e n and x be positive): 



^(x, £ri) 



exp (—e n x/u 4- x/2l ee ) 



x 2 Fi[l/2 + e„, 1/2, 1; x(x)} , (3.17) 



where 2Fi(a, b, c; z) is the hypergeometric function, 

X (x) = 1 - exp(2x/; ee ), £„ = £ n /2vT . (3.18) 

We now analyze the asymptotic behavior of &(x,e n ) 
as a function of two dimensionless parameters x/l ee and 
£ n /2TrT. For x/lee > 1 and e n /2TrT > 1, Eq. $3l7\i gives 



exp ( — £ n x/u - x/2l e e) ^(^n) 



iwir \/uv 



(3.19) 



where T(z) is the gamma- function. After the analytical 
continuation to real energies i£ n — > e + iO, Eq. (|3.19|) 
reveals oscillations of G\(x,e) as a function of £x/u 
and an exponential decay as a function of x/l ee - Using 
Eq. (|3.19[) for the analytical continuation is only accurate 
for (\e\/T)(x/l ee ) 3> 1. In the opposite limit, one has to 
analytically continue already in Eq. (|3.17p . which yields 
the "static limit" for the Green's function with 



2exp(x/2/ ee ) r -, 
1^'°) - -—7=r. K [x{x)\ 



VK\/UV 



(3.20) 
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where K(z) is the complete elliptic integral. For x/l ee 3> 
1, Eq. reduces to 



<f(x> Z ee ,0) 



inJuV l ee 



■ exp(—x/2l e 



(3.21) 



Continued to real energies, Eqs. (|3.19[) and (|3.21[) match 
onto each other at (\e\/T)(x/l ee ) ~ 1. Finally, the high- 
energy short-distance asymptotic behavior of the Green's 
function is given by 

9{* « U*. » T) * eXp( - £ ^ /M+) I ( ^ 



(3.22) 

where Jo (2) is the Bessel function of the imaginary argu- 
ment. 

The imaginary part of G r + (x,e) [obtained as an ana- 
lytical continuation of G\(x, is n ) onto the real axis of 
energy from the upper half-plane] as a function of x for 
small and large e/T is shown in Figs. [5] and respec- 
tively. While in the former case there are only simple os- 
cillations which are suppressed exponentially on the scale 
of l ee [cf. Eq. (|3.19p ], the behavior of G r + [x, e) in the lat- 
ter case is richer. Specifically, Fig. [6] exhibits beatings 
and an intermediate power-law decay, in agreement with 
Eq. (|3.22p . The real part of G r + (x, e) behaves similarly. 



uv ImG+(x, e) 




FIG. 5: Imaginary part of the Green's function G r + {x,e) [ob- 
tained as an analytical continuation of Eq. (|3,16[1 onto the 
real axis of e from above] as a function of x for small energies 
£ < T shows oscillations cos(xe/u) exp(—x/2l ee ) with a pe- 
riod 2-ku/e, exponentially suppressed on the scale of l ee - The 
parameters of the plot are: e/T — 0.25 and u/v = 1.1. 



C. Green's function in the (q,e) representation: 
Spectral weight 

Here, we complete the analysis of the single-particle 
Green's function in a spinful LL by inspecting its spec- 
tral properties in the momentum-energy representation. 
Fourier-transforming Eq. (|3.14[) with respect to x, and 
analytically continuing the result onto the real axis of 
e from the upper half-plane, ie n — > e + iO, we get the 
retarded Green's function G+(q,e) in the (q, e) represen- 
tation. Similarly, the analytical continuation onto the 



uv ImG + (x, e) 
1 




FIG. 6: Imaginary part of the Green's function G r + (x,e) as 
a function of x for large energies e S> T. One sees oscilla- 
tions with a period 2nu + /e, beatings with a period 2nu-/e, 
and a power-law decay a; -1 / 2 (taking place up to x ~ l ee ). 
The asymptotic behavior for large x S> l ee at which G r + (x, e) 
is suppressed exponentially is not shown in the figure. The 
parameters of the plot are: e/T = 20 and u/v — 1.1. 



real axis from below, ie n — > e — zO, yields the advanced 
Green's function G+{q, e). 

The retarded and advanced Green's functions of right 
movers can be written in the form 



2L 



V{±K U )V{±H V ) 



(3.23) 



where the signs + and — correspond to the retarded (R) 
and advanced (A) Green's functions, respectively, 



(e/u—q)l ee: K v = (e/v — q)l e 



and 



V{z) = 



r[(l - 2zz)/4] 
r[(3-2w)/4] 



(3.24) 



(3.25) 



As a function of complex variable z, V(z) has a series of 
simple poles (originating from the gamma-function in the 
numerator) at z = — i(4m + l)/2, where m is a positive 
integer. The pole that is closest to the real axis [corre- 
sponding to complex e = uq—iu/2l ee and e = vq—iv /2l ee 
in Eq. (|3.23p ] determines the spatial/temporal decay of 
the Green's functions considered in Sees. IIII Al and UlI Bl 
An alternative representation, which straightforwardly 
splits G^j A {q,e) into the real and imaginary parts, is 



Gf A (q,e) 



I, 



-R,A 



where 

C R ' A (x,y) = smh[n(x + y)/2] 
=F zcosh[7r(x — y)/2] 

and the real function K(z) is given by 
1 



K(z) 



2tt 



T[(l-2iz)/A]T[(l+2iz)/A] 



(3.26) 



(3.27) 



(3.28) 



The upper and lower signs in Eqs. (|3.26|) — (|3.28|) corre- 
spond to the retarded and advanced functions, respec- 
tively. 

Note that the influence of interaction on the behavior 
of G+' A (q 7 e) is twofold. Firstly, it factorizes the single- 
particle fermionic Green's function into two parts charac- 
terized by different velocities — v for the spin factor and 
u for the charge factor. This is the essence of the SCS. 
Secondly, at finite T, interaction leads to a broadening 
of the singularities in the spectral weight, i.e., to a shift 
of the singularities of G+' A (q,e) into the complex plane. 
In the (x,e) representation, this shift manifests itself in 
the exponential damping of the Green's function on the 
spatial scale of l ee , as discussed in Sec. MI 61 

Figures [7J and [5] illustrate how the imaginary and real 
parts of G+(q,e) as a function of e evolve with varying 
temperature. At T = one gets^ 



G*(q, e) = (3.29) 

■s/e-uq^s-vq 

(for u — > v, e is understood as e + iO). At low T, there 
is a double-peak structure which represents the SCS, 
with square-root singularities at e — vq and e = uq, 
weakly smeared by temperature. With increasing T, 
the broadening becomes more pronounced and eventu- 
ally two peaks in the spectral weight merge into a single 
peak of width ~ aT . 

At this point, it is worth recalling that we have ne- 
glected effects of interaction which are related to the 
exponent ctb ~ 0(a 2 ) in Eq. (|3.2p . Retaining a^, i.e., 
including the last factor in Eq. (|3.2[) would only lead 
to the following two effects, both of which are of mi- 
nor importance in our consideration at weak interaction. 
Firstly, there will be an additional small asymmetry be- 
tween two peaks in Figs. [7] and [51 Secondly, an additional, 
weak singularity (characterized by the exponent cvb) will 
arise at e = — uq (cf. Refs. HH and HI, where the single- 
particle Green's function in a spinful LL was investigated 
at T = beyond the weak-interaction limit). 

To conclude, in Sec. IIIII we have analyzed the behav- 
ior of single-particle excitations in a spinful LL at finite 
T. We have demonstrated that, because of the SCS, it is 
dramatically modified as compared to the spinless case. 
In particular, the decay length l ee turned out to be para- 
metrically shorter than for spinless electrons. However, 
a priori it is not immediately clear to what extent the 
modification of the single-particle properties will affect 
the transport (i.e., two-particle for fcrmio ns) pro perties 
of a spinful LL. Indeed, as shown in Refs. [27ll28l for the 
spinless case, the WL dephasing length l^, is parametri- 
cally longer than l ee . Calculation of the conductivity of 
a disordered quantum wire in the presence of spin is a 
subject of the next section. 
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-Im G*(q,e) 




vq uq 



FIG. 7: Spectral weight ImG+(</, e) (arbitrary units) for right 
movers as a function of energy for different temperatures and 
u/v — 1.1. At low T, the SCS manifests itself in the double- 
peak structure, with two peaks located at e = vq and e = uq. 
When T increases (T/vq = 0, 0.05, 0.1, 0.5 and 1) the singu- 
larities are rounded off and eventually a single peak emerges. 

-Re G*(q,e) 





I 






vq ^ 


r 



FIG. 8: Real part of the right-mover Green's function 
G+(g, e) as a function of energy. Parameters are the same 
as in Fig. [7J The singularities are smoothened and the peak- 
dip structure broadened with increasing temperature. 

IV. WEAK LOCALIZATION 

So far we have analyzed the single-particle properties 
of a spinful LL at finite T in the absence of disorder. 
Now we introduce disorder and turn to the calculation 
of a two-particle quantity, namely, the conductivity. At 
high T, the leading term in the conductivity is given by 
the Drude formula, 

2e 2 

o-u = 1, (4.1) 

with a renormalize d 56 ' 57 i 58 ' 59 by interaction — therefore 
temperature dependent — mean free path I oc T Qr . 

We consider now a correction Act to the Drude con- 
ductivity, associated with the quantum interference of 
electron waves multiple-scattered off disorder. In ID, 
the leading contribution to Act comes from a Cooperon- 
type scattering process which, in contrast to higher di- 
mensionalities, involves a minimal possible number of 
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(a) 




(b) (c) 




FIG. 9: Diagrams giving the leading contribution to the in- 
terference correction to the conductivity. The dashed lines 
represent the impurity scatterings and the solid lines denote 
the electron Green's functions (with the disorder effects incor- 
porated at the self-energy level) and the dashed lines repre- 
sent the impurity-induced backscattering. The diagrams are 
understood as "dressed" by interaction as shown in Fig. \W\ 



scatterings on impurities, namely three ( "three- impurity 



Cooperon" ). 27 i 28 The peculiarity of ID in this respect is 
that a single-channel quantum wire is in the WL regime — 
not strongly localized — only if the dephasing length 
that cuts off the WL correction is shorter than the mean 
free path I. That is, the WL correction is accumulated 
on ballistic scales — hence the shortest possible Cooperon 
ladder with three impurity legs. 



A. General expression for Cooperon 

The leading term in Act is given by the diagram o 27 ! 28 
in Fig. [51 These are understood as dressed by interaction- 
induced fluctuating gauge factors exp[±iO ^{x , t)] at- 
tached pairwise to the backscattering vertices, as de- 
scribed in Sec. Ill Cl and illustrated in Fig. [10] for the case 
of diagram (a). The sum of contributions of diagrams 
(b) and (c) is equal to that of diagram (a). At this level, 
there is no difference in the structure of the diagrams be- 
tween the spinful and spinless cases — the only difference 
stems from the particular form of the correlators of the 
phases 9^{x, t). 

Averaging over the fields O^x, r), we get for the in- 
terefernce correction at Matsubara frequency O m : 



T 



Acr(iO TO ) = 2 x 2 x 2 x (evf x ( ^— ) x -!- 

l/T ,1/T ,1/T rl/T ,1/T ,1/T f f f 

dri / dfx \ dr 2 / df 2 / dr 3 / df 3 / dx\ \ dx 2 / dx 3 
to Jo Jo Jo Jo Jo J J J 

x [g+(xi - £3,7-1 - f 3 )Q' 1 (xi - X3,n - f 3 )] [g-(x 2 - xi,t 2 - t 1 )Q~ 1 (x 2 - x\,t 2 - n)] 

x [g+(x 3 - x 2 ,t 3 - t 2 )Q~ 1 (x 3 - x 2 ,t 3 - t 2 )] [g-(xi - x 3 ,n - t 3 )Q' 1 (xx - x 3 ,?i - t 3 )] 

x [g+(x 2 - xi,f 2 - f 1 )Q~ 1 {x 2 - xi,f 2 - n)] [g-(x 3 - x 2 ,f 3 - f 2 )Q~ 1 (x 3 - x 2 ,f 3 - f 2 )] 

x Q(x! - x 3 ,n - t 3 ) Q(xi - x 3 ,fi - f 3 ) Q(xi - X 2l Tl - f 2 ) 

x Q(x 2 - X\,T 2 - n) Q(x 3 - X 2 ,T 3 - f 2 ) Q(x 2 - X 3 ,T 2 - f 3 ) 

x Q-\0, n - n) Q- x (0, t 2 - f 2 ) Q-^O, t 3 - f 3 ) 

x Wlixi -x 3 ,T U f 3 ,Q m )yvl(xi -xa,n,T 3 ,Cl m ) , (4.2) 

where L is the system size, the free Green's functions g±(x, t) are given by Eq. (|3.ip . the interaction-induced factors 
Q(x, t) are defined by Eq. (|2.22p . and the factors W % ± (x, r, r', J7 m ) come from integration of the two Green's functions 
attached to the current vertices over the external coordinates and times: 



WUx,t, r',Cl m ) 



sgnf2 r , 



;{(« 



-j:o m T 



-iVt m T ! 



O m | + v/l 
e- in ™ T 9(n m x) 



) + \n m \i ( 



-\Q m x\/v y 



(4.3) 



r 



w!: f (x,T,T',n m ) = -w l + f (~x,T,T',n m ), 

Wl(x,T,T',Sl m ) = W l + (x,T,T',-n m ). 



(4.4) 



In Eq. ()4.3j) . we have included vertex corrections for the 
current vertices, which arise from the anisotropy of im- 
purity scattering [recall that only backscattering off im- 
purities (|2.3[) is retained in the model]. The vertex cor- 
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in the limit of small a is then written as 




FIG. 10: The same diagram as in Fig. [9^, with the interaction- 
induced factors exp(±jS,j) shown explicitly. The solid lines 
with arrows stand for bare Green's functions, the crosses 
for the impurity vertices, and the wavy lines for the factors 
exp(±20 M ). The space-time coordinates of the backscattering 
vertices are denoted by N = (xn,tn) and N = (xn,tn)- 
Averaging over the fields 9^ couples all wavy lines with each 
other, cf. Fig. H 



rections result in a replacement of the total scattering 
rate v/2l by the transport scattering rate v/l. Note also 
that Eq. (|4.2[) is written in terms of the contribution of 
the Cooperon with chiralities of the current vertices as 
shown in Fig. [TUJ The numerical coefficient in Eq. (|4.2|) 
takes into account that Acjwl is a factor of 2x2x2 = 8 
larger than the contribution of the diagram in FiglTO] [one 
of the factors of 2 comes from the spin, another from a 
summation over chiralities of the current vertices and the 
third one from diagrams (b) and (c) in Fig. [9] . 

The approximation (|3.4j) means that the terms in the 
exponent of Q(x,t) [Eqs. (j2~2"T|) . d2~22]) . and([2J26])] that 

come from -B++ (x, r) and B (x, r) and are proportional 

to a b are neglected. As for the term that comes from 

B_i {x,t) and is proportional to a r , it leads to a renor- 

malization of the impurity strength (see Ref. |28 for de- 
tails) but does not contribute to the dephasing rate to 
first order in a, similarly to the spinless case. Therefore, 
in the calculation below we put both a b and a r in Q(x, t) 
equal to zero, while Iq is replaced by I which is understood 
as the renormalized mean free path. The factor Q(x,t) 



Q(x,t) 



< sinh[7rT(:r/u + it)] sinh[7rT(:r/u — it)] 
sinh[7rT(a;/i; + it)] sinh[7rT(a;/w — it)] 



(4.5) 



The second term in Eq. (|4.3[) . proportional to I /I, can 
be omitted for l^/l <C 1, so that one more approximation 
we make is to take the factors (|4.3[) in Eq. (|4.2[) at x = 0: 



\n m \ + v/i 



It is convenient to introduce new variables 



(4.6) 



and 



3 (1 


= Xl 


- x 3 




= x 3 


- x 2 


x c 


= Xi 


- x 2 


T a 


= n 


- f 3 


T c 


= r 2 


- n 


n 


= T 3 


- r 2) 




= n 


- t 3 




= r 2 


- n 


n 


= r 3 


- f 2 



(4.7) 



(4.8) 



These satisfy the constraints x c — x a + x\> and r a + r;, + 
t c + f a + fb + f c — 0. We thus represent Eq. (|4.2p as an 
integral over the variables (|4.7p , (|4.8p and insert, instead 
of T/L, in the integrand the factor 



-Tl, + T,. + T„+T h + T, 



] 8{x c ~x a -x b ) (4.9) 



which contains summation over fermionic Matsubara fre- 
quencies. By extracting the vertex functions (|4.3|) at 
x = 0, Eq. (|4.2p in the limit a <§; 1 is rewritten in the 
new variables as 



~ lim \ - — — — - y~] / dT a dT b dT c df a df b df c / dx a dx b dx c 

x G + (x a ,T a ) G-(x a ,f a ) G + (x c ,t c ) G-(x c ,f c ) G + (x b ,n) G-(x b ,fb) 

x C + (x a ,T a ) C-(x a ,f a ) C + (x c ,t c ) C-(x c ,f c ) C + (x b ,n) C-(x b ,fb) 

X Q(x a ,Th +T c )Q{x ai f b +fc)Q(x c ,T a +f b )Q(x c ,T a +T b )Q(x b ,f a +f c )Q(x b ,T a + T c ) 

x exp{i(tt m +e n )(T a + T b + T c ) -ie n (? a +T b + T c )] 5(x a + x b - x c )\ , (4-10) 
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where G±(x, r) is given by Eq. (|3.4p . Q(x, r) by Eq. (|4.5p . and we introduce 



C±(x,t) = 



sinh[7rT(ai/i> =F : 



sinh[7rT(a;/u =p it)] 

such that g±Q~ 1 — G±C±. When deriving Eq. (14. 10[) we have used the approximation (14.61) 



w;(o,T 1 ,f 3 ,fi m )w{(o,f 1 ,r 3 ,r!, 



i 



(\n m \ + v /i) 2 

_|_ e if2 m (Tl-T 3 ) I e «n m (f3-Tl) 



in m (ri — ti) _ „ifi m (T3— T3) 



2 gjn m (T 3 -T 3 ) 



(4.11) 



(4.12) 



The "diagonal" terms exp[iO m (fi — ti)] and exp[if2 m (T 3 — 
f 3 )] yield identical contributions [which is accounted for 
by the factor of 2 in Eq. (|4. 1 0[l ] ■ The cross-terms are 
neglected, since, after the integration over times (see 
Sees. lIVBl and IVl below), they produce the products of 
Green's functions in the (x, e) representation of the type 

G^ a) (x,e n )G r _ (a \x,e k )=0 (4.13) 

[see Eq. (I3.15[) ]. This corresponds to retaining only those 
Cooperon diagrams that contain an equal number of the 
retarded and advanced Green's functions, even with e-e 
interaction included. 

We are thus left with the functions G±(x,t), Q(x,t), 
and C±(x,t) that are all given by various combina- 
tions of square-root factors y / 'smh[TrT(x/u ± it)] and 
y/ smb.[irT(x/v ± it)]. Note that G±(x,t) and C±(x,t) 
enter Eq. (|4.10p only in the combination G±C± with the 
same arguments. Equation (|4.10[) will be analyzed in the 
next two subsections. 

B. Regular impurity configurations: Weak 
localization 

In Eq. (|4.10[) , we transform the integration contours 
for each of the time variables similarly to Fig. |H As 
a result, we obtain integrals along square-root branch 
cuts in the vertical direction, each of which connects two 
points whose coordinates can be written as t + ix/u and 
t + ix/v with different r and x. For example, let us 
assume, here and throughout the paper below, that in 
Eq. ((1710]) all 

x a ,x b ,x c >0 (4-14) 

(the region of integration x a < and Xb, x c > gives the 
same contribution to Actwl)- Then, starting with the 
integration over r a and closing the contour upwards, we 
represent the integral along the real axis of r a as a sum of 
three integrals around vertical cuts: between ix a /u and 
ix a /v, between —fj, + ix c /u and — fj, + ix c /v, and be- 
tween — f c + ixb/u and — r c + ixb/v. The first cut comes 



from the Green's function G + (x ai T a ), whereas the last 
two from the factors Q{xb, T a + r c ) and Q{x c , T a + 77,), re- 
spectively. Since G+{x,t) as a function of x for a given 
r falls of on the scale of It, while Q(x, r) on the scale of 
l ee , one sees that in the limit of small a the main con- 
tribution to Actwl comes from the branch cuts that are 
associated with the Green's functions. The cuts related 
to the factors Q(x,t) can be neglected. 

The selection of singularities at a <C 1 is closely anal- 
ogous to that in the spinless case, see Appendix F of 
Ref. For spinless electrons, the main contribution to 
Actwl stems from singularities ( "nearly poles" for weak 
interaction) of the single-particle Green's functions at the 
classical trajectory of an electron moving with the veloc- 
ity u. Other close-to-pole singularities (those in the fac- 
tors Q and C), which are related to "nonclassical trajec- 
tories" , yield subleading corrections small in the parame- 
ter It /l<j> <C 1. The spinful problem is very much similar 
in this respect. The main difference is that the dominant 
singularities are now pairs of close square-root branching 
points rather than the poles. The transformation of the 
poles into the branch cuts, induced by the SCS, can be 
viewed as a "smearing" of the classical trajectories: all 
velocities between v and u become accessible. 

Let us first consider the contribution to Eq. (|4.10p from 
typical impurity configuration for which the characteris- 
tic scale of 

x a ~ x b ~ x c (4-15) 

is of the order of It- We can then expand all sum's 
in the C and Q factors as 

sinh(TrTy) ~ ie 7 ^^ 1 sgny , \y\ > 1 . (4.16) 

It is immediately seen that using Eq. (14. 16)) reduces the 
product of six C factors to a simple exponential: 

C ~ exp(2a; c /; ee ) . (4.17) 

Similarly, the product of six Q factors 
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Q = Q 
x Q 

is represented as 
where 



x a ,i \ t c + t b + 

Xc: 1 [ ~t~ 



2x b ~\~ X a 

u 

X b x a 

u 



X a ,i \ t c + t b + 

Xb,i [t a +t c + 



*2xfr -|- x a 
u 

*2x a -\- Xfr 
u 



Q 



3C6, M ia + *c + 



X a X b 

u 

^X a -\- X b 

u 



Q ~ exp(7rTg/2) 



(4.18) 
(4.19) 



(j = \2xb/u + t c + tb\ - \2xb/u + t c + tb — 2x a /u-\ + \2x c /u + t c + tb\ - \2x c /u + t c + h + 2x a /u_| 

+ \2xb/u + t c + tb\ — \2xb/u + t c + t b — 2x a /u-\ + \2x c /u + t c + ib\ — \2x c /u + i c + tb + 2x a /u_| 

+ \2xb/u -t a + ib\ - \2xb/u - t a +tb + 2x c /u-\ + \2x a /u + t a - t b \ - \2x a /u + t a — tb + 2x c /u-\ 

+ \2x a /u + t a - t b \ - \2x a /u + i a -tb + 2x c /u-\ + \2x b /u - i a + t b \ - \2x b /u- i a + t b + 2x c /u_| 

+ \2x a /u + t a + t c \ - \2x a /u + t a +t c ~ 2x b /u-\ + \2x c /u + t a + t c \ - \2x c /u + t a + t c + 2x fc /u_| 

+ \2x a /u + t a + t c \ - \2x a /u + t a + i c - 2x b /u-\ + \2x c /u + i a + t c \ - \2x c /u + i a + i c + 2x b /u^\ . (4.20) 

In Eqs. (|4.18p and (|4. 20[) we have shifted the time variables according to 

ITj = -tj-Xj/u, 



iTj = tj +Xj/U, 



(4.21) 



with j = a, b, c. The shifted variables tj and tj in Eq. (|4.21[) are real on the branch cuts corresponding to the Green's 



functions in Eq. (|4. 1 0[) and change along these cuts from to 2xj/i 



j/u, where u~ is given by Eq. (|3.1ip . 



Inspecting Eq. (14.201) . we observe that for x a > 2x c u/u- (or similarly for Xb) all the moduli can in fact be omitted 
on the Green's function branch cuts, which yields 

Q ~ exp(-47rTx c /u_) = cxp{-2x c /l ee ) (4.22) 

(the opposite case x a ^ < 2ax c is addressed in Sec. fV|. As a result, the factors C and Q compensate each other: 

QC ~ 1 . (4.23) 

The integrand of Eq. (|4.10p thus reduces to a product of the single-particle Green's functions in the (x, e) representation 
[Eqs. ;E3. (I3~T71) ]: 



A(T WL r 2ttT 

= lim < — 

ctd n-»o I 



/>oo 

>J / dx a dx b dx c 5 (x a + x b - x c ) 
„ Jo 



n m (\n m \ + v/i) 2 

x G\_ (x a , ie n + if2 m )GI!_ (x c , is n 4- iVL m )G\(xb, is n + ifi m )G™ (x a , ie n )G" (x c , ie n )G" (x b , is n ) \ , 

(4.24) 

where we have taken into account that only terms with e n < 0, e n + Q m > survive, in view of Eq. (|4.13|) . which 
follows from Eq. J3.15[) . Performing the analytical continuation to real frequencies f2 m — * fl + iO, we get in the dc 
limit fl -> 0: 



CD 



-A 



WL 



T 



where the numerical factor v4wl is defined by 

dz 



Awl = 7T 



cosh ttz 



dx 



dyK{x, z)K(y, z)K(x + y + xy, z), 



with 



K(x, z) = 2 Fj ( 1/2 + iz, 1/2, 1; -x ) 2 Fj ( 1/2 - iz, 1/2, 1; -x ) 



(4.25) 

(4.26) 
(4.27) 



We have estimated ^4wl — 0.13 by taking the integral (|4.26p numerically. 



The small factor (l ee /l) 2 in Eq. (|4.25|) is due to the exponential decay exp(— 2x c /l ee ) of the product of six 
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Green's functions in the integrand of Eq. (|4. 24[) on the 
scale of l ee . One sees that the dephasing factor that sup- 
presses the interference term in the conductivity behaves 
as exp(— Lc/l$), where Lq = x a + Xb + x c — 2x c is the 
total length of the Cooperon loop and the WL dephasing 
length 1$ reads 

h = lee ■ (4.28) 
Schematically, Eq. (|4.24[) can be estimated as 



(a) 



J ^Y-exp(-2x a /l ee ) 

f°° dXb flee\ 2 

J — cxp(-2x b /l ee ) ~ -y-j-J 



(4.29) 

so that for typical realizations of disorder with x a ~ x b ~ 
x c each of the distances is of the order of l ee , see Fig.fTTh. 
For comparison, in the spinless case ; 27 ' 28 the relevant dis- 
tances obey x a Xb ~ ll ee , which yields 



Arr sl 
a(7 WL 

0"D 



f dx a f dx b 
- / — j— / —;- exp{-x a x b /ll ee ) 
Jo 1 Jo ' 



lee n I 



(4.30) 

It follows that for spinful electrons is much shorter 
than the dephasing length for spinless electrons (recall 
that for the latter, diverge o 27 i 28 in the limit of vanish- 
ing disorder). In fact, in the spinful case 1$ is equal to 
the single-particle (electron) decay length l ee , in contrast 
to the spinless case, where 1$ » l ee . 



V. ANOMALOUS IMPURITY 
CONFIGURATIONS: MEMORY EFFECTS 

The WL contribution to the conductivity, calculated 
in the preceding section, is associated with scattering on 
compact three-impurity configurations which are "regu- 
lar" in the sense that the characteristic distances between 
all three impurities are the same. Below, we will see 
that "anomalous" (strongly asymmetric) configurations 
in which two of the impurities are anomalously close to 
each other, i.e., x a <§; x b or Xb <C x a (see Fig. [Tib), 
give rise to a conductivity correction which is larger than 
Actwl if T is sufficiently high. As mentioned already in 
Sec. HI the relevance of the asymmetric configurations is 
related to the classical ME eo in electron kinetics, in con- 
trast to the quantum interference of scattered waves that 
yields Actvvl- 



A. Qualitative consideration: Identifying scales 
and parameters 

To demonstrate the peculiarity of the asymmetric im- 
purity configurations, it is instructive to consider first 



(b) 



FIG. 11: Three- impurity configurations that give the main 
contribution to the (a) quantum (WL) and (b) classical (ME) 
corrections to the conductivity. The characteristic distances 
between the impurities are shown. 



the limit of two scattering events occurring at the same 
point, by setting 



X\ = x 3 



(5.1) 



in Eq. (|4.2| . As discussed at the beginning of Sec. II VB . 
for typical impurity configurations (|4.15p the main con- 
tribution to the correction (|4. 2[) comes from "smeared" 
classical trajectories, meaning that all trajectories with 
velocities between v and u contribute to Actwl- The 
case (I5.ip is, however, special in that only one velocity 
remains and that is the velocity of nonintcracting elec- 
trons v. Indeed, at x\ = X3 the four Q factors that de- 
pend on x\ — X3 drop out of the integration over times — 
since Q(0, r) [Eq. (|4.5[) ] does not depend on r. After 
this, the integrals in Eq. (|4. 2|) over t\ and f\ are dom- 
inated by the poles of the noninteracting Green's func- 
tions g+(0, n — T3) and <?-(0, f\ — T3), which yields 



ti = r 3 , 



ri 



73, 



(5.2) 



and then all the remaining factors Q compensate each 
other. We thus end up with a product of six bare Green's 
functions g± multiplied by the factors VV± , which consti- 
tutes the noninteracting limit of the problem. 61 

In fact, the compensation of the dephasing factors 
is evident even before the averaging over the fluctu- 
ating fields 6±(x,t) (Fig. ITOj) . Indeed, the factors 
exp{±i[6>+(a;, t) — 9-(x, r)]}, dressing two impurity ver- 
tices, cancel each other when taken at the same space- 
time point [Eqs. (|5.ip and (|5.2p ]. As a result, the third 
impurity located at x-i becomes decoupled with respect 
to e-e interaction from the two-impurity complex at 
x\ = x 3 . 

In the cyclic variables (|4.7p . Eq. (|5.ip means x a = 0. 
Since interaction completely drops out of the problem at 
x a = 0, the integration over the remaining spatial coordi- 
nate x b — x c in Eq. (|4.10p is not cut off by dephasing, in 
contrast to the regular impurity configurations, for which 
it is restricted to x c < l ee <C I- In this situation, we 
have to take into account a disorder-induced damping 
of the single-particle Green's functions g±(x, r), which 
results in an additional factor exp(— |a;|/4Z) attached to 
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FIG. 12: Each diagram for Act in Fig. |9]has two ME contri- 
butions, coming from x a ~ It (3 — > 1, 1 — > 3) and Xb ~ Zt 
(2 — f 3, 2 — ► 3). The remote-impurity line should be replaced 
by the full diffuson ladder. 



each g±(x,r). On the Cooperon loop, these combine 
to give the overall factor exp(— x c /l) in the integrand 
of Eq. (H]), so that then limited by the mean 

free path, see Fig. [TTb . 

Since at x a — the characteristic distance to the re- 



mote third impurity in the three-impurity Cooperon di- 
agram happens to be in the crossover region between the 
ballistic and diffusive motion, wc should extend the single 
scattering on the third impurity to an infinite sequence 
of scatterings on other impurities i.e., to a diffuson 
ladder, as shown in Fig. [T21 One sees that the diagram 
takes the form characteristic of a quasiclassical an 
electron is scattered at x\ ~ X3, then moves around dif- 
fusively, and returns to X\ ~ £3 where it is scattered once 
again. Clearly, this is a non-Markovian process which is 
beyond the conventional Boltzmann description. How- 
ever, the non-Boltzmann type of kinetics associated with 
the return processes is classical in origin: as discussed 
above, when the points x\ and X3 are sufficiently close to 
each other, dephasing becomes irrelevant. We will first 
analyze the simplest three-impurity diagram at finite but 
small x a and include the diffusive returns (which will only 
rcnormalize the numerical prefactor) later in the end of 
this section. 

It is instructive to begin with a semi-quantitative anal- 
ysis which will give correctly the parametric dependence 
of the result but not the numerical coefficient. To this 
end, we replace all hyperbolic sines in both the Green's 
functions G± and the Q and C± factors in Eq. (|4.10|) by 
their exponential asymptotics, Eq. (|4.16[) . This approx- 
imation is parametrically correct, since all integrals in 
Eq. (|4.16|) are determined by the regions of integration 
in which the arguments of the hyperbolic sines are of the 
order of or larger than 1. Within the "exponential" ap- 
proximation, the product C of six factors C'± is given by 
Eq. (|4.17[) . while the product Q of six Green's functions 
in Eq. (|4. 16|) becomes 

Q cx exp(-2x c /l ee ) , (5.3) 

so that the exponential factors in Q and C cancel each 
other, similarly to the regular impurity configurations 
in Sec. HVBl What is different, however, is that when 
calculating the factor Q [Eqs. (|4.18[1 - (|4.20[) ] one can no 
longer omit the moduli in Eq. (I4.20|) if 

x a < 2ax c (5-4) 

[or x < 2ax c — we will proceed with the estimate for the 
case of Eq. (|5.4p ]. For the strongly asymmetric configu- 
rations (|5.4[) , Eq. (|4.20[) is rewritten as 



q ~ — 8x c /u^ + \2x a /u — t \ — \2x a /u — t + 2x c /ii- 1 + \2x a /u — if, | — \2x a /u — t + 2x c /u- 1 

+ \2x a /u + t c \ — \2x a /u + t c — 2x c /u- J + \2x a /u + i c \ — \2x a /u + i c — 2x c /u- 1 , (5-5) 

where we have neglected terms of second order in a. In particular, the integration over t a and i a goes along very 
short branch cuts of length oc ax a ~ a 2 x c . To our accuracy, the cuts reduce to poles. For this reason, we have set 
t a = la = in Eq. J53J). 

The integrals over remaining time variables in Eq. (|4.10p decouple from each other and can be easily calculated 
(note that the integrals over t and t c are identical — in the dc limit — to those over ib and i c , respectively). Let us 
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denote Acr asym the contribution to the conductivity that comes from the strongly asymmetric impurity configurations 
with \x a b\ < 2ax c . Carrying out the analytical continuation for both fermionic and bosonic frequencies and taking 
the dc limit, we get, within the exponential approximation: 



Act, 



0"D 



de dx a dx c 
Tcosh 2 ( e /2T) 

9(x a — ax c )0(2ax c — x a 



exp ( — y> exp 



2x t 

I P P 



2lx 2Z ee ttT 



X a ) 



I X a X c c 

h 1 2/? 2/7' 



e 

7TT 



X a X c £ \ 



2Z T ' 27 



(5.6) 



where the functions /i,2,3 are defined as follows 

cxp(izy) — 1 



f 2 (x,y,z) = exp[iz(y- 
exp(izx) 



. exp(«za:) — 1 



i z 



— i exp(a; — y) 



exp[(iz + l)(y — x)] — 1 



fa(x,y,z) 



1 



+ z exp(— x) 



{iz + l) 

exp[(iz + l)y] — exp[(i.z + l)a;] 



(5.7) 



(5.8) 



{iz + l) 

In Eq. (|5.6p . the 9- functions split the domain of integration over x a into two, ax c < x a < 2ax c and 

< x a < ax c . 

In the first interval of x a , the integration over x c is limited by the factor exp(— 2x c /l ee ) in the first line of Eq. 
which yields a contribution to ACT aS y m / ctq of the order of —l ee lT/l 2 - This is much smaller than the contribution of 
the regular impurity configurations [Eq. (|4.22j) ] and can be neglected. However, the situation is qualitatively different 
for the interval (|5.8[) . Indeed, for l,x^y, the functions \ f 2 {x,y, z)\ 2 and \ f 3 (x,y, z)\ 2 take the form 



\f 2 (x,y,z)\ 2 
\f 3 (x,y,z)\ 2 



exp(izx) — 1 



1 



exp[2(y-x)] 



(5.9) 
(5.10) 



As a result, the factor exp(— 2x c /l ee ) from the first line of Eq. (|5.6j) is canceled by the factor exp(2y) from Eq. (|5.10|) . 
The remaining integral over x c is no longer restricted to x c < l ee but rather extends to x c ~ I, in agreement with our 
qualitative consideration of the ME effect for x a — at the beginning of this section. The ME contribution Actme 
to the conductivity thus arises from the impurity configurations obeying Eq. (|5.8j) and is given by the second term in 
Eq. (|5.6p . It can be estimated as 



Act me 

CD 



f°°dx c ( x c \ [ ax <dx a ( 2x a 

Jo — exp l--H — exp rv 



It 
' I 



and becomes much larger than the WL contribution [Eq. (|4.22p ] if T is sufficiently high, namely if 

T-»v/a 2 l . 



(5.11) 



(5.12) 



The characteristic regions of x a ,x that control the quantum (WL) and classical (ME) corrections to the conductivity 
are shown in Fig. 113a . For comparison, an analogous plot for the spinless problem is presented in Fig. 113b . One 
sees that in the spinful case the "quantum" and "classical" domains are parametrically separated from each other. 
As a result, their areas l 2 e oc T~ 2 and Ut oc T -1 , respectively, can "compete" with each other. On the contrary, in 
the spinless case, the ME strips of area ~ ll ee are directly adjacent to the WL domain, whose area ~ ll ee \n(l/l ee ) 
is (logarithmically) larger. It follows that for spinless electrons the ME correction Actj^ e gives only a subleading 
contribution as compared to the WL correction Act|^ l [Eq. (|4.30| ]: 



Act 8 ' 



<7D 



Act s ' 



0"D 



In 



(5.13) 



[recall that 1$ ~ (llee) 1 ^ 2 and l ee ~ a 2 It for spinless electrons]. 
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(a) (b) 




FIG. 13: Domains of the distances x a ,xt that govern the WL (gray-shaded area) and ME (dashed area) contributions to the 
conductivity for the spinful (a) and spinless (b) models. 



B. Rigorous calculation of Acme 



In the preceding subsection, we have adopted the "exponential approximation" by replacing all the hyperbolic sines 
in the correlation functions by their asymptotics (|4.16p . This allowed us to identify the relevant scales and obtain the 
parametric dependence of the result. Here, we calculating the integrals over the temporal variables tb, tb, t c , and t c 
more accurately and find the numerical prefactor in Eq. (|5.11|) . The estimate in Sec. IV Al teaches us that the dominant 
contribution to Actme comes from x a ~ It and x b ~ x c ~ I, while the integrals over tb lC and ib >c are determined by 
the upper limit of integration, i.e., by the time scale ax c /v ~ al/v. In the WL regime, this time scale is much larger 
than 1 /T. 

Similarly to the case of regular impurity configurations, we observe that the product of all the C factors can always 
be approximated by Eq. (|4.17[) , whereas the hyperbolic sines should be retained in the Green's functions G±. As 
for the factors Q, in contrast to the regular configurations, not all of them can be replaced by their asymptotics. 
Specifically, those hyperbolic sines in the Q factors that correspond to the moduli in Eq. (|5.5[) (i.e., all terms in q 
except for the first one) should retain their form and not be replaced by the exponentials. Importantly, within this 
"partially exponential" approximation, all integrals over times are still decoupled. Thus, we replace the exponential 
factor depending on t c in Sec. IV Al according to 



](x a /l T + nTt c - x c /l ee ) - i8( - x a /l T ~ nTtc + x c /l ee ) 



x exp 



- {-x c /l ee + \x a /l T + irTt c \ - \x a /l T + irTt c - x c /l ee \) - |e„ + fl m \ t c 



2exp(-|e„ + f2 m | t c 



sinh 1 ^ 2 (a; a //T + TtTt c ) 



sinh 1 ^ 2 (7rTi c ) smh 1 ^ 2 {x c /l ee — nTt c ) sinh 1 ^ 2 (a; a /^ + ^Tt c — x c /l ee ) 
2exp(x a /2l T ) exp(-|e„ + f2 m | t c ) 



sinh 1/2 (a; c // ee - irTt c ) sinh 1/2 (a; a /Z T + irTt c - x c /l ee ) 
and make a similar replacement for t c (with s n + Q m — > e n ). The factors depending on tb are modified as follows: 

1 



(5.14) 



i exp 



{-x c /l ee -I- \x a /l T - TrTt b \ - \x a /l T - nTt b + x c /l ee \) - \e n + Q m \ t b 



2exp(-\e n + n m \t b ) 



sinh 1 ^ 2 (xa/lr — TtTtb) 



sinh 1/2 (7rTt b ) sinh 1/2 (a; c /; ee - itTt h ) sinh 1/2 (x a // T - itTt b + x c /l ee ) 
2exp (—x a /2l T ) exp {-\e n + O m | tb) 



sinh 1/2 (a; c // ee - nTt b ) smh 1/z (x a /l T - nTt b + x c /l ee ) 



1/2/ 



(5.15) 
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the integral over i& again differs only in that e n + VL m changes to e„. Using Eqs. (]5 . 14[) and (|5. 15|) . we get 

= -A MB — , (5.16) 

(td t 

where ^4me is given by the dimensionless integral (check the first arguments of F) 

f°° dz f°° ,2 ,2 

j4me = tt / 2 — / dx \M(z,x)\ \M(z,—x)\ . 



cosh ttz Jo 

(5.17) 

The function M(z,x) in Eq. (f57T7|) is defined by 

1 f^, exu(2izy) 



1/2 (y) sinh 1/2 (y + a;) 



We have estimated ^4me numerically as Ame — 0.2. 



Finally, let us discuss the overall combinatorial factor 
in Actme (which is already included in ^4me)- Firstly, 
similarly to the WL correction, the contribution of dia- 
gram (a) in Fig. [12] should be multiplied by a factor of 
2x2x2 = 8 due to (i) two possible chiralities of the 
current vertices; (ii) diagrams (b) and (c) in Fig. [T2l and 
(iii) two possible anomalous configurations for each di- 
agram: x a ~ It <C X(, and Xb ~ It "C x a . Secondly, 
as mentioned above, not only the return after one single 
backscattering but rather the entire diffuson ladder con- 
tributes to Actme [Fig. H2] • The insertion of the diffuson 
into the three-impurity diagram effectively generates an 
additional velocity-vertex correction, which yields a fac- 
tor of 1/2 (the ratio of the transport and total scattering 
times for the backscattering impurities). 



VI. PATH-INTEGRAL METHOD 

So far, we have been treating the conductivity cor- 
rection within the formalism of Matsubara functional 
bosonization. This powerful method treats on an equal 
footing 2 ^ the real inelastic scattering processes responsi- 
ble for dephasing and the virtual transitions responsible 
for the renormalization effects. An alternative approach, 
formulated for the spinless case in Refs. [2?] and l28l con- 
sists of two steps. Firstly, disorder is renormalized by vir- 
tual processes with characteristic energy transfer larger 
than T. What is obtained after the renormalization is 
an effective "low-energy" theory which is free of ultravi- 
olet singularities characteristic of a LL. The low-energy 
theory is treated by means of a path-integral approach, 
analogous to the one developed in Ref. [63| for higher- 
dimensional systems. This method is particularly con- 
venient for an analysis of inelastic scattering (dephas- 
ing) in problems with a nontrivial infrared behavior. In 
Ref. HI, the WL correction to the conductivity was cal- 
culated both within the path-integral and the functional- 
bosonization schemes, with identical results. The path- 



I 

integral calculation also allows one to "visualize" the ori- 
gin of the dephasing processes in terms of quasiclassical 
trajectories. 

In this section, we present a path-integral analysis of 
the conductivity correction in the spinful problem. It 
turns out that the situation here is more intricate than 
in the spinless case, for two reasons: (i) the character- 
istic energy transfer in the path-integral calculation is 
of the order of T (while it was much smaller than T for 
spinless electrons); (ii) because of the SCS, the velocity of 
the quasiclassical trajectories is not uniquely defined: the 
whole interval of velocities between v and u contributes 
to Actvvl, see Sec. IIII1 As a result, we will only be able 
to reproduce the parametrical dependence of the conduc- 
tivity correction but not the numerical prefactor. Never- 
theless, this analysis is useful, since it yields a physically 
transparent picture of the quantum interference and de- 
phasing in the spin-charge separated system. 

Since the dominant contribution to the dephasing rate 
is given by the (^-processes of scattering between elec- 
trons from the same chiral branch, we will set g% = 
which is consistent with our central approximation ctb — > 
0. Furthermore, it is convenient to introduce different 
coupling constants ji* 2 - g\ and g±, for interaction between 
electrons with equal and opposite spins, respectively (at 
the end, we set g\ = g± = g ~ 2irva). In view of the 

Pauli principle, the g\ -interaction does not lead to any 
real scattering but only renormalizes the velocity of elec- 
trons: 

v -> v* = v(l+g\/2irv) . (6.1) 

This allows us to put g\ — 0, simultaneously replacing v 
by v*. All the nontrivial interaction- induced physics is 
due to g^- -processes. 

Solving the corresponding RPA equations for the in- 
teraction propagators, we get for parallel (V^) and anti- 
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backward 




forward 




forward 



FIG. 14: World lines corresponding to the time-reversed tra- 
jectories (solid and dashed lines) with velocity v* in the three- 
impurity Cooperon. The dash-dotted line describes the prop- 
agation of a plasmon with velocity u. (a) Regular configu- 
ration (Sec. IIV|) contributing to Actwl; (b) asymmetric con- 
figuration (Sec. |V| contributing to Actme- The interaction 
line gives a contribution to the dephasing action S which is 
proportional to a(Nf — Nt) 2 . Here Nf t t is the number of the 
small-angle intersections (black dots) of a plasmon line with 
the forward (/) and backward (6) paths. The intersections at 
a large angle (the plasmon and electrons moving in opposite 
directions, unfilled circles) contribute to the dephasing only 
at second order in a and hence are neglected. One sees that 
Nf 7^ Nb for typical configurations (a), whereas for asymmet- 
ric configurations (b) Nf = Nb for most of the plasmon lines. 
For a corresponding plot in the spinless case, see Fig. 7 of 
Ref . HI 



parallel (V ) spins : 



V^ + (q,Lo) = 2irv*a 



(v*q — u)Y 



(uq — u>)(yq — lu) 



Vl + (q,u) = -2irv 



*„2 



v*q(v*q — u>) 
(uq — u))(vq — lo) 



and 



(6.2) 



(6.3) 



where u = v(l + 2a) and v* = v(l + a). We note that 
V does not enter the path integral, since electron spin 
is conserved. 

The general expression for the WL dephasing ac- 
tion acquired on the exactly time-reversed trajectories 
x b(t) — X{(tc — t) for the three- impurity Cooperon 
reads 2 ^ 



Suite, {x(t)})-. 



-T 



to r-tc 
dti / dt2 

n Jo 



~2H 



dq 
2^ 



x 



\mV} v {q,u) 



exp{iq[xi(ti) - x 3 (t 2 )} - iu[h - t 2 ]} 

UJ 

= -T I C dh f ° dti^Mh) - Xj (t 2 ), ti-ta] , 
o Jo 



(6.4) 

where each of the indices i,j takes one of the values "f" 
(for the forward path of the Cooperon) or "b" (for the 

backward path), \x — sgni^ and v — sgnij, and J 7 ^, (x, t) 

is the Fourier transform of uj^^mV^uiq^)- The main 
contribution comes from the diagonal terms with i = j 
and t\ = t%, for which fj, = v, see Fig. [14k . The imagi- 
nary part of the corresponding interaction propagator is 
written as 



ImVl ± (q,uj) = {irv*a)' 



to 



x [ uS(vq =F uj) + v 5(uq =p w)] , (6.5) 



which gives 



F± ± {x,t) 



{v*af 
2uv 



[uS(x =fvt)+v S(x =F ut) } 



(6.6) 



A graphic illustration of the e-e scattering processes 
contributing to the dephasing action S = 2(Sg — Sfb) 
is presented in Fig. 1141 There, we show by the circles 
and black dots the space-time coordinates for which the 
arguments of the (^-functions in Eq. (|6.6[) are zero. The 
main contribution to S comes from small-angle intersec- 
tions (black dots). For electron trajectories characterized 
by velocity v* , these give a large factor of the order of 
a -1 , either v*/\u-v*\ or v*/\v-v*\ [cf. Eq. (|2~14|) ]. As 
a result, the action for typical impurity configurations 
(Fig. [T3k) reads 



S - a 2 Tt c /\a\ ~ \a\Tt c 



(6.7) 



For comparison, in the spinless cas o 27 i 28 , intersections 
between the interaction and particle propagators run- 
ning in the same direction give zero dephasing because of 
the cancellation between the Hartree and exchange terms 
(Sec. Ill Bj) . On the other hand, at large angles the ballis- 
tic (straight line) interaction propagator always intersects 
a pair of (forward and backward) electron trajectories, 
which gives zero dephasing as well. As a result, only the 
intersections at large angles that arise from scattering of 
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the interaction propagator off disorder contribute to the 
dephasing rate. 27 ' 28 

More rigorously, substituting Eq. (|6.6p in Eq. (|6.4| , we 
find the dephasing action for the trajectory characterized 
by the velocity v*, 

!x a , < x a < ax c /2 , 

ax c /2, ax c /2 < x a < (1 — a/2) x c , 
x c -x a , (1 - a/2) x c < x a < x c , 

which is illustrated in Fig. 1151 Any other velocity be- 
tween v and u yields a qualitatively similar action, but 
with different numerical factors. This is another indica- 
tion of the previously discussed fact that all velocities 
from this interval contribute to the conductivity correc- 
tion. 

s 



ME WL ME 




a/2 l-a/2 1 



FIG. 15: Dephasing action S (in units of 2x c /lee) correspond- 
ing to the velocity v* , Eq. (|6.8|l . as a function of x a for fixed 
x c . 

One sees from Eq. (|6.8p (middle line) and Fig. [T5] (flat 
region) that for typical impurity configurations the de- 
phasing action becomes of the order of unity at tc = 
2x c /v* ~ 1/aT. This gives a parametric estimate for 
the dephasing length 1$ ~ a^ 1 ^, and for the WL cor- 
rection Actwl/cd ~ ~{lee/l) 2 i m agreement with the re- 
sults of Sec. lIVl Equation (|6.8[) (first and third lines) and 
Fig. [15] also demonstrate the suppression of the dephas- 
ing action in the anomalous (strongly asymmetric) im- 
purity configurations (Fig. [Tib) responsible for the ME, 
Acme/cd ~ —It/I, as discussed in Sec. El 

VII. SUMMARY 

To conclude, we have analyzed, within the functional 
bosonization formalism, the quantum interference of in- 
teracting electrons in a disordered spinful LL. Our results 
are summarized in Table [Tj and in Fig. 1161 

The single-particle properties of fermionic excitations 
in this model have been studied in Sec. IIIII in several 
representations. Two most important and interrelated 



features of the single-particle spectral characteristics are: 
(i) the SCS, as a result of which the whole range of veloc- 
ities between the charge velocity u and the spin velocity 
v contributes to the spectral function, and (ii) the single- 
particle decay on the spatial scale of l ee , Eq. (|3.9[) . 

In Sec lIVI we have calculated the leading quantum in- 
terference correction to the conductivity, Eq. (|4.25j) . The 
corresponding dephasing length 1$ is given by Eq. (14.28)) . 
Two qualitative differences as compared to the spinless 
case should be emphasized in this context. Firstly, for 
spinful electrons, the decay length of single-particle exci- 
tations is inversely proportional to the first order in the 
interaction strength a, l ee ~ a~ l lT, whereas l ee ~ a~ 2 lx 
in the absence of spin. Secondly, the WL dephasing 
length 1$ is equal to the single-particle dephasing length 
l ee in the spinful model, whereas 1$ depends on the 
strength of disorder, ~ \Jl ee U f° r spinless electrons. 

In Sec. [V] we have analyzed the contribution to the 
Cooperon diagrams of "nontypical" configurations of dis- 
order with two impurities being anomalously close to 
each other. We have shown that this contribution, 
Eq. (|5.16[) . describes the quasiclassical ME. It is para- 
metrically smaller than the WL term at sufficiently low 
T and gives the leading correction to the conductivity in 
the limit of high T. 




vll vial v/a 2 l 

FIG. 16: Schematic plot (log-log scale) of the T dependence 
of the conductivity er(T) of spinful electrons. 

Our findings, demonstrating the strong dependence of 
the quantum interference effects on spin in a LL, imply 
that the Zeeman splitting by magnetic field should lead 
to strong effects in the conductivity of a single-channel 
quantum wire. Results obtained in this direction will be 
reported elsewhere^ 2 - 
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useful discussions and gratefully acknowledge support by 
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ter for Functional Nanostructures, and by RFBR [Grant 
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e-e scattering 
length l ee 


dephasing 
length Ifi 


WL correction 
Actwl/ctd 


ME correction 

A<7me/<Td 


spinless 


V 




(Uy, I vln(a 2 Tl/v) 
U ) U l4, a 2 Tl 


lee v 

I a 2 Tl 


spinful 


V 

aT 


V 

aT 


-&)'--(£,)' 


It v 
I Tl 



TABLE I: Characteristic spatial scales induced by interaction and the conductivity corrections for spinless and spinful disordered 
LLs. 
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